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Resume 

Nous elaborons une methode originale d'imagerie directe d'un objet. II s'agit d'une 
methode de type derivee topologique, destinee a inverser des mesures contenant 
l'amplitude des ondes retrodiffusees a la suite d'illuminations laser. 

Abstract 

Laser imaging. We derive an original direct imaging method of an object. It is 
based on topological derivatives, and aims at inverting the amplitude of waves that 
are retropropagated after laser illuminations. 



Introduction 



L' identification d'objets par un systeme d'imagerie laser est un sujet aux do- 
maines d'application varies, comme la detection d'engins explosifs improvises 
(IED) ou de tumeurs dans le derme. Dans ce contexte motivant, nous elabo- 
rons une methode d'imagerie directe (non iterative) par un tel systeme. Ainsi, 
nous developpons un algorithme d'inversion de mesures d'ondes retrodiffu- 
sees par l'objet a identifier, a la suite d'illuminations laser. Cette technique 
est monofrequentielle : une seule frequence est exploitee ; elle est monosta- 
tique : les sources d'illumination et les recepteurs sont confondus. Comme 
en pratique la vue est partielle (une partie de l'objet n'est pas illuminee), 
nous supposerons que les mesures sont prises sur une droite. Enfin, comme 
pour les applications, nous supposerons que les mesures contiennent l'ampli- 
tude (ou l'intensite), mais ne contiennent pas la phase. Comme les methodes 
d'inversion directes usuelles (MUSIC [I], backpropagation method, Kirchhoff 
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migration) reposent sur la connaissance de la phase, le principal interet et 
l'originalite de notre methode concerne ce point. 

Pour ajouter de l'information, de fagon non restrictive par rapport aux appli- 
cations, nous supposons que l'on peut prendre les mesures sur une autre droite, 
parallele a la premiere. Ainsi, pour chaque angle d'eclairage, nous disposons de 
deux mesures d'amplitude prises dans la meme direction. En approximation 
haute frequence, le passage d'une mesure d'amplitude a l'autre s'obtient par 
une equation de transport ; ceci cache une certaine notion de phase dans le 
rapport entre les deux mesures. L'idee est de trouver les petits disques, qui, 
eclaires, retrodiffusent une onde dont la phase permet de « transporter » les 
premieres mesures vers les deuxiemes. En plus d'utiliser directement les me- 
sures, cette idee permet d'exploiter l'information cachee au sujet de la phase. 
Cette methode peut etre formulee a l'aide d'un cout que l'on cherche a rendre 
petit. La « validite » d'un point est alors teste par la valeur de la derivee to- 
pologique du cout : un point est bien place par rapport a l'objet recherche 
si cette derivee est tres negative. II s'agit ainsi d'une approche par analyse 
asymptotique topologique, dans le meme esprit que celle conduite dans [3f2] . 



1 Position du probleme 

1.1 Probleme direct. On considere l'espace libre M 2 avec permittivite 
6q > et permeabilite fj, > 0. On ajoute dans ce milieu une inclusion D de 
permittivite > et permeabilite yU* > 0. D est supposee pres de l'origine. 
On eclaire D par un laser, i.e., une onde incidente plane d'angle unitaire t], 
de pulsation u > : u- mc (r],x) = e lkov ' x . Le nombre d'onde associe est ko := 
Uy/eoJIo, et la longueur d'onde du laser est A := Le nombre d'onde associe 
a D est fc* := uo^Je*n*. L'onde qui se propage est alors u[D](rj,x) = u(x), 
solution de : 



avec une condition de radiation sortante sur l'onde diffractee u^a [D](r], x) : = 
u[D] (?7, x) — u mc (r], x), x G M 2 \ D. Numeriquement, ce probleme est resolu par 
la methode des equations integrales, comme suggere dans 

1.2 Probleme inverse. On fixe deux droites de mesures paralleles en 
champ lointain : T : x 2 = 7o e ^ Ti : ^2 = 7i, °h l es constantes 70, 71 > sont 
grandes par rapport a la longueur d'onde, avec 70 — 71 raisonnable (ni trop 
grand, ni trop petit). En eclairant l'objet D avec l'angle x/\x\, l'amplitude 
retrodiffusee en x est A[D](x) = |Mdif[-D](x/ \x\ , x)\, notee A [-D] et ^i[-D] 
respect ivement sur r et Ti. On suppose qu'un objet D* est situe a proximite 
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de l'origine, avec pour parametres electromagnetiques > et > connus. 
On suppose que l'on mesure A [-D*] et AJ-D*]. Le probleme inverse que l'on 
souhaite resoudre consiste a retrouver l'objet a partir des mesures, en 
inversant A[-D*] par A, i.e. en resolvant A[D] = A[D*\, sur T U Fx. 



2 Methode de derivee topologique 

2.1 Demarche. Par approximation haute frequence, on approche l'am- 
plitude Ai[JD#] en transportant l'amplitude Ao[-D«], avec une certaine phase 
y[.D*] : si A satisfait l'equation de transport : V A ■ Vip[D*\ + ^AA^D*] = 
0, avec A = Aq[D*\ sur r , alors la trace de A sur r\ approche Ai [£)*]. 
D'autre part, en transportant A [-D*], avec la phase connue (p[D e (z)] issue 
de l'onde diffractee par un disque D e (z) de centre z et de petit rayon e > : 
VA ■ Vip[D E (z)] + ±AA(p[D £ (z)} = 0, avec A = A [D*\ sur T , alors on ob- 
tient sur r x une amplitude t4i[D £ (^)]. On a alors l'intuition que z est un 
point bien positionne par rapport a l'objet lorsque l'amplitude transpor- 
ts t4i[D £ (z)] est proche de l'amplitude mesuree A\ [D*\. L'idee originale de 
reconstruction que nous proposons consiste a chercher les z qui rendent petite 
la quantite ||Al[_D £ (z)] — Ai[-D*]||. Nous obtenons une expression asympto- 
tique de AJ.D e (2:)] par la methode des caracteristiques. Puis un developpe- 
ment asymptotique de || AJ.D e (,z)] — Ai[-D*]|| permet d'identifier une derivee 
topologique a rendre tres negative pour estimer D*. 

2.2 Derivee topologique. La resolution de l'equation de transport par 
la methode des caracteristiques conduit au theoreme suivant. 

Theoreme 2.1 Pour tout p £ r , en posant t] = I 'angle d'eclairage, le 
transport de la mesure A [D*](p) vers T 1; le long de la direction de retropro- 
pagation (z,r]) ; avec phase (p[D E (z)](ri, ■), conduit a l'amplitude : 

MD.w M = M»,M (^f J" 2 (i - ' ^.^.-J + ° ^ 

oil x v = z + 71 ~ Z2 1] = Ti fl (z, rj), et t^ = — 71 ~ 70 est le temps de trajet. 

Parametrons par x = (xi,7i) : r](x) = ^E^r\> p{ x ) = z + ^rff ( x ~ z ) > e ^ P osons 

a {x,z) := Ao[D*M*)) i^f'^ a ^ z ) : = -\ a ^ Pour un 

domaine D, notons Al[.D] l'amplitude issue du transport de Ao[-D*] sur Ti, 
dans la direction de retropropagation, avec la phase ip[D] de l'onde diffractee 
par D. Le theoreme precedent justifie le resultat asymptotique suivant : 

Theoreme 2.2 (Derivee topologique) La derivee topologique de I'applica- 
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tion D i — > \\Ai[D] — Ai[D*]\\, evaluee en z, dans la direction le disque unite, 
est : 

d(z) = rr— r f , .. (a (- , z) - A X [D*\, z)) . 

\\a (-,z) - Ai[D^\\\ 

2.3 Methode d'inversion directe. Si le petit disque D £ (z) est bien 
place, i.e. « pres » de alors la phase tp[D £ (z)] approche la phase V't-D*]? 
et done v4i[D e (z)] approche Ai [£)*]. L'approche originale que nous proposons 
pour reconstruire la forme consiste a minimiser D i— > \\Ai[D] — Ai[D*]\\, 
parmi les petits disques. Le developpement asymptotique du theoreme : 

\\A x [D e {z)] - A 1 [D,)\\ = |M-, z) - + ed(z) + O (e 2 ) , 

suggere que ||Ai[£) e (z)] — est petit lorsque la derivee topologique d(z) 

est negative et tres grande en valeur absolue. Pour estimer D*, on trace ainsi les 
lignes de niveau de z i — > d(z), puis on en choisit une de niveau do arbitraire. 
Celle-ci peut servir de point depart pour un algorithme iteratif. 
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